Noise-induced phase transitions in field-dependent relaxational dynamics: Tlie 

Gaussian ansatz 
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We present an analytic mean field theory for relaxational dynamics in spatially extended systems 
that undergo purely noise-induced phase transitions to ordered states. The theory augments the 
usual mean field approach with a Gaussian ansatz that yields quantitatively accurate results for 
strong coupling. We obtain analytic results not only for steady state mean fields and distribution 
widths, but also for the dynamical approach to a steady state or to collective oscillatory behaviors in 
multi-field systems. Because the theory yields dynamical information, it can also predict the initial- 
condition-dependent final state (disordered state, steady or oscillatory ordered state) in multistable 
arrays. 
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I. INTRODUCTION 



The interplay of stochasticity and nonlinearity often 
leads to a wide range of interesting and often counter- 
intuitive behaviors. Noise and nonlinearities have been 
implicated as sources of both zero-dimensional transi- 
tions [3| as well as nonequilibrium phase transitions in 
spatially extended systems [2,]. Particularly intriguing 
among these is a broad class of macroscopic ordering 
phenomena that have been called "purely-noise-induced" 
because order is exhibited only in the presence of noise. 
These include order-disorder transitions [1, 0, Q, the 
emergence of periodic spatial structures and pattern for- 
mation in increasingly noisy systems [1, 0], and, more 
recently, purely noise-induced collective oscillatory be- 
havior [Si, J3, lOj . Purely noise-induced phase transitions 
to oscillatory behavior may be dynamically induced and 
dependent on the Stratonovich drift, i.e., on the inter- 
pretation of the noise 3,, ^] , or it may occur in systems 
with field dependent relaxational dynamics where the in- 
terpretation of the noise at most shifts the transition 
points 0, d, [13] ■ The particular novelty of the latter 
transition lies in the fact that it arises entirely from an en- 
ergy functional-like relaxational dynamics, and not from 
short-time dynamic instabilities which become strength- 
ened and sustained by spatial coupling. As such, these 
models represent systems that, while decidedly nonequi- 
librium and potentially time-dependent in nature, nev- 
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ertheless maintain certain parallels with equilibrium sta- 
tistical models, namely the evolution towards the mini- 
mum of a (possibly time-dependent) effective potential. 
Our focus in this paper is on phase transitions in systems 
with field dependent relaxational dynamics. 

Analytic work on these models has been fairly lim- 
ited, with most of the information coming from numer- 
ical simulations. We are not aware of solvable analytic 
theories for locally coupled (e.g. nearest neighbor) noisy 
extended systems in which any sort of ordering transition 
is observed as the system or noise parameters are var- 
ied. Globally coupled arrays have been more amenable 
to mean field theories based on the analytic solution of 
a Fokker-Planck equation in the stationary state. This 
approach has been successfully applied to systems that 
achieve a time-independent steady state. Such mean 
field theories provide information about the nature of the 
steady state and the conditions that lead to disordered 
vs ordered states. Thus, while quite useful for ascertain- 
ing asymptotic properties, these static theories provide 
no insight into the dynamical evolution toward a steady 
state. Such insight is particularly important in multi- 
stable systems. 

Most strikingly, these theories fail to provide even a 
reasonable asymptotic description of phase transitions to 
collective time-dependent behavior, when the complexity 
of the corresponding Fokker-Planck equation prohibits a 
self-consistent analytical treatment even in globally 
coupled arrays. The complicated behavior seen, for ex- 
ample, in the transition to noise-induced limit cycles thus 
requires different methods for analytic study |^]. 

Our purpose here is to provide an approximate analyt- 
ical method supplementary to the mean field approach 
for studying the dynamics of noise-induced phase transi- 
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tions in relaxational systems inspired by a related quan- 
tum field method ^lli] . By introducing a Gaussian ansatz 
and a series expansion about the (time-dependent) mean 
field values, we develop a time-dependent theory which 
captures the dynamics of these systems in the limit of 
strong spatial coupling. As well as providing new ana- 
lytical insight into the evolution to the steady state in 
one-field systems, this treatment yields a set of simple 
approximate ordinary differential equations detailing the 
oscillatory dynamics and the approach to these dynam- 
ics in two-field systems. We note that an alternative ap- 
proach to the one adopted in this paper is to generalize 
the procedure developed in Q that gives a set of differen- 
tial equations for the central moments. It can be shown 
that this alternative approach gives the same results at 
leading order for large coupling as the Gaussian ansatz, 
we are preparing a manuscript where we show the details 
of the comparison. 

The paper is organized as follows. In Sec.|lT]we briefly 
recall the one- and two-field relaxational models. The 
mean field equations and the Gaussian ansatz for their 
solution are presented in Sec. IIIII In Sec. IIVI we analyze 
the one-field system and test our Gaussian ansatz results 
in the steady state and in the dynamical approach to 
the steady state, including regimes of multistability. The 
two-field system is examined in Sec. |Vl where we show 
that our theory captures the collective oscillatory behav- 
ior, including regimes of multistability. In both Sees. IIVI 
and |V] we explore approximate solutions to the differ- 
ential equations obtained with the Gaussian ansatz to 
successfully reduce the problems to quadrature in some 
regimes. We conclude with a short summary in Sec. I VII 



II. THE MODELS 

We briefly recall the one- field and two-field relaxational 
models. In the one-field case we introduced the Langevin 
equation defined on a lattice 



[n^,)Y'^at). (1) 



where Lpi is the value of the scalar field at site i, 
is the field-dependent kinetic coefhcient, J-{{tf}) is an 
energy functional, and £,i{t) is a spatio-temporal white 
noise with zero mean and intensity cr^, {£,i{t)^j{t')) = 
a-^SijS{t — t'). All the quantities introduced previously 
and those that we introduce later in this paper are con- 
sistently adimensionalized. As the choice of noise in- 
terpretation does not qualitatively affect the dynam- 
ics d, H, we choose the Ito interpretation for con- 
venience. The energy functional includes local poten- 
tials Vjifii) and a simple harmonic coupling between 
sites 
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Here A'' is the number of lattice sites and K is the cou- 
pling strength. The sum runs over all n sites j 
coupled to site i. For nearest neighbor coupling n = 2d 
while for global coupling n = A^ — 1 (other forms of cou- 
pling were also considered @ but will not be dealt with 
here). In [sj we carried out comprehensive studies of 
the phase space diagram in {K, ct^) space where different 
disordered and ordered phases can be found for different 
types of local potentials V{(p) and of kinetic coefficients 
r((p). In general, necessary conditions on these functions 
for the observation of order-disorder phase transitions in- 
clude that, at increasingly large \(p\, V{(f) grow at least 
harmonically, and T{(p) go to zero so that the effect of the 
noise is weakened. We were able to analyze the steady 
states for these systems analytically within a mean field 
approximation. This approximation leads to a Fokker- 
Planck equation that can be solved in the steady state, 
which allowed us to establish the phase space diagram for 
the steady state, but provided no dynamical information. 

The addition of a second field to this problem was in- 
troduced in fio'l, 

^, = ^Tif,,)^-^M1 + [r((p,)]^/'e.(0 - coz,, 

ofi (3) 

where a; is a frequency. This system undergoes a noise- 
induced phase transition to collective oscillatory behav- 
ior when the noise exceeds a critical intensity. The mean 
field Fokker-Planck equation for this two-field system can 
not be solved analytically because it remains time depen- 
dent for all time. Thus, all the information obtained up 
to now has been numerical. 

Since our Gaussian ansatz method relies on the mean 
field evolution equations, in the next section we briefly 
review the mean field approach and introduce the ansatz. 



III. THE MEAN FIELD AND THE GAUSSIAN 
ANSATZ 

In the mean field approximation the sum over neigh- 
bors connected to site i in the derivative of the energy 
functional appearing in Eq. ^ or Eq. ([31), 



K 



is replaced by the mean field value. 



(4) 
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Since all the sites are then equivalent, the lattice index 
can be dropped and the set of field equations reduces to 
a single equation (one- field system). 



^(t) =a(^;(^o(<)) + [r((^)]^C(i), 



(6) 
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or to two coupled equations (two-field system), 



the Langevin equation 



Z = LOip, 



where 



(7) 



(8) 



The unknown mean field lySolO niust be determined self- 
consistently, 



(9) 



Here 



stands for the statistical average with respect 



to the probability density p associated with Eq. ([6]) or 
Eq. 0. 

The Fokker-Planck equation for the probability den- 
sity ^((p, t\ (po{t)) in the one- field case follows immediately 
from the Langevin equation ([6|). 
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We explicitly note the dependence of p on the unknown 
mean field </'o(0: which must be determined via Eq. ([9]) 
using the solution of the Fokker-Planck equation. The 
time dependent probability density has not been found 
analytically. In the stationary state the left side of 
Eq. ([To]) is set to zero and the equation can be solved, to 
yield the steady state probability density 



Pst(¥';<^o) = ^(<^o)r((^) 

= iV(^o)e-('/'^')^=^^(^) 
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where 



y^sAv) ^ y^^) + f (^0 - + Yinr(^) (12) 

and N(po) is the normalization constant. The mean field 
(po can then be found from Eq. ([9|) . A disordered station- 
ary phase is associated with the solution (po = 0, while 
a solution po ^ corresponds to an ordered stationary 
phase. The phase boundaries for different forms of V{(p) 
and T{p) are detailed in f^. In particular, it shows that 
even if the potential V{p) is monostable (single- well), 
the presence of field-dependent multiplicative noise leads 
to an effective bistable (double-well) potential in certain 
regions of the parameter space. This is precisely the man- 
ifestation of a noise-induced phase transition. More elab- 
orate forms of r((^) could even lead to multistable poten- 
tials. Note that the procedure in [5] leads to a complete 
mean field stationary state analysis, but does not provide 
information about the dynamics of the approach to the 
steady state. 

For the two-field case the Fokker-Planck equation for 
the probability density p(p, z, t; po{t), zoit)) follows from 
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This equation has not been solved analytically. Further- 
more, since we know from numerical simulations that 
the system supports collective oscillations, it is necessary 
to solve the time-dependent problem even to find the 
long-time behavior. The only available information to 
date is numerical. 

To obtain time-dependent solutions to the one-field 
and two-field models, we will assume a Gaussian form 
for the evolving probability density with time-dependent 
parameters to be found self-consistently from the asso- 
ciated Fokker-Planck equation. Thus, for the one- field 
problem we take 



p{p,t) 



(14) 



where the time dependent mean field (po(0 ^'^^ inverse 
width parameter A(t) are to be determined. The param- 
eter C{t) is found from the normalization condition 



dpp{p) = 1 



to be given by 



C{t) 
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-A{t) 



(15) 
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It is immediately evident from a comparison of Eqs. (jl4|) 
and pT|) that is at best an approximate solution, 
but we will subsequently show that for sufhciently strong 
coupling K this in fact provides an excellent approxima- 
tion for the mean value and width of the distribution. 
Note that Eq. pT|) implies that in the mean field approx- 
imation A = — Vgjj((/?o)/c2 for the stationary state. For 
the two-field case we posit the form 



p{p,Z,t) — gMv-Va)'^ + E{v-Va)(z-za)+M{z~-ZQf+C^ j-^^y-j 

where now 29(^)7 E{t) and M{t) are also to be deter- 
mined. The normalization condition 



fixes C(t), 



dpdzp{p, z, t) = 1 



^ 1 fAMA-E^ 
^=2^H 



(18) 



(19) 



Note that the Gaussian only has finite norm for 4M A — 
>0. 

We expect these forms to work best for strong cou- 
pling. If coupling is too weak, then there is no transition 
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to collective behavior and a Gaussian ansatz is not ap- 
propriate. As coupling becomes extremely strong, the 
distribution approaches a i5-function. This provides the 
motivation for a narrow distribution, whose mean and 
width we assume to be well captured by a Gaussian when 
coupling is strong. The Gaussian ansatz thus rests on the 
observation that ipo{t) ~ zo{t) ~ A{{S(p)'^) ~ E{S(pSz) ^ 
M{{5zY) ~ 0{K^) together with the assumption that 
the second order moments {{5(pY) ^ {{^z)^) ^ 0{1/K) 
and the expectation that higher order moments are sub- 
leading for large K. Here Sip = ip — ipo and Sz = z — zq. 
In the following sections we show that these assumptions 
lead to consistent results. 



IV. SINGLE FIELD 

We start with Eqs. (fT4|) and (fTB]) . substitute the Gaus- 
sian into the Fokker-Planck equation pil)l , Taylor expand 
a, r, and their derivatives, and implement a large- if ex- 
pansion considering the previous comments. We also rec- 
ognize that the kinetic coefficient and its derivatives, as 
well as the local potential and its derivatives, are inde- 
pendent of K. The contributions to leading orders in K 
then result in the set of equations 



A ^ 2ATo{K + (7^A) 
^0 = 



- ( ^ + f ) n, 



(20) 
(21) 



For each function f{ip) we have adopted the notation 
/o = f{^o), fo = [df{>f)/dip]^^^„, etc. Note that it is 
clear from these two equations that the system evolves 
towards A ~ 0{K) and (po ~ 0(1). We also point 
out that the normalization condition (jl6p is consistent 
with the evolution equation one obtains for C{t), namely, 
C = To{K + a'^A). The set of equations ^ and ^ is 
of course nonlinear and can not be solved exactly analyt- 
ically, but it is merely a set of two ordinary differential 
equations whose numerical solution is trivial. 

Even so, we can make further analytic progress by not- 
ing that for large K, and provided the initial value of A is 
of (negative) 0{K), the evolution of A toward its station- 
ary state is clearly faster than that of (/jq, thus allowing us 
to consider Fq as nearly constant during the relaxation of 
A. In this approximation Eq. (|20[) is a Riccati equation 
with solution 



A{t)^ 
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After the stationary state for A is reached, the evolution 
of po{t) is governed by 



d(po 



(23) 



Before comparing the results of the theory with those 
of numerical simulations, we note that in the stationary 
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FIG. 1: (Color online) The time evolution of ifio (top panels) 
and A (bottom panels), that give the mean and the inverse 
width of the distribution, is shown for simulations of glob- 
ally coupled arrays of N — 4000 sites and compared to the 
numerical solution of Eqs. (|20p and (I2ip . The left column 
represents modest coupling, K — 10, while the right column 
shows K = 80; = 3.5 for all plots. The initial values 
~A{0) are chosen to be > K/a^; specifically, A{0) = —9.7 in 
the left column and A{0) — —58 in the right column. The 
dark (blue) lines show data from lattice simulations, while 
the light (red) lines represent theoretical predictions. The 
very light (brown) curves in the lower right panel and inset 
represent the uncoupled dynamics given by Eq. (|23p and the 
Riccati equation (|22p . The lower left inset shows the sim- 
ulation results and the prediction of the mean field theory 
(solid black line) in the stationary state, which is exact when 
N —> oo. The inset in the bottom right figure shows a close 
up of the late time evolution. Note the difi^erent horizontal 
and vertical scales in the various panels. 



state where A — ip = we find from Eq. or Eq. 
that A = —K/a'^, and that ipo obtained from Eq. pT|) or 
is the solution of the condition dVeff{'Po)/dipo — 0. 
The latter is exactly the mean field solution to the prob- 
lem, which is thus recovered from the Gaussian ansatz. 
The former differs from the exact inverse width of the 
distribution PT|) by contributions of 0(1). 

However, the ansatz takes us beyond the stationary 
solution to provide information about the dynamics of 
the system as it approaches the steady state. In Fig. [T] we 
show four sets of results for the mean field foit) and the 
inverse width parameter A{t). One is the outcome of the 
direct simulation of globally coupled arrays for moderate 
and for strong coupling K. The second is the result of 
integrating Eqs. ([^0)1 and (PT|) . the third is the outcome 
of Eqs. and that assume different relaxation 
rates, and the fourth is the outcome of the mean field 
distribution pT|) . In these and all subsequent figures we 
have made the representative choices 



V{ip) = 



r(^) 



i + ip^ 



2' ''''' l + ip^' 
which were also used in our earlier work [^, [lo| . 



(24) 



The first conclusion is that the time scale of relaxation 
and the steady state value of (pQ are correctly predicted 
by the large coupling theory, even for modest values of 
K. The agreement is spectacular during the approach 
to steady state. In addition, we note that the time scale 
separation between A and (po is quite evident for large K 
and is a reasonable assumption also for modest values, 
K ^ 10. Further discussion of the results in this fig- 
ure requires that we take note of the different horizontal 
and vertical scales in the different panels. A small but 
consistent discrepancy between the Gaussian theory re- 
sults for A and simulation results for the inverse width of 
the distribution arises at long times owing to neglected 
higher order contributions in the theory. In the lower 
left panel we have expanded the vertical scale to make 
the difference clear, but note that it is extremely small 
for large K (see lower right inset). The lower left inset 
confirms that the simulation and mean field theory sta- 
tionary widths are in fact identical, as they should be 
for sufhciently large arrays. While we only show results 
for two values of K and one lattice size, the results for 
a range of values of K {K = 20, 40, 60), of lattice sizes 
(N = 250 up to = 16000), and of initial conditions 
in which A{t = 0) —K follow the patterns described 
above. We thus conclude that the Gaussian ansatz theory 
gives quantitatively excellent results for the evolution of 
the globally coupled system toward the steady state and 
for the steady state itself when the coupling is strong and 
the width of the initial distribution is of the same order 
as that of the steady state. 

Our theory is based on a mean field theory, and so the 
appropriate comparison with simulations is as we have 
shown in Fig. [2 with a globally coupled array. How- 
ever, mean field theories are often used to describe lo- 
cally coupled systems, and so we compare our theoreti- 
cal results with simulation results in which the units in a 
two-dimensional array are connected only to their near- 
est neighbors. The results are shown in Fig. [51 The 
theory accurately captures the behavior of ipo, includ- 
ing both the transient and steady state dynamics, in the 
large and modest K regimes. However, the inverse width 
of the distribution (as given indirectly by A) is under- 
estimated by the theory, as it is by the original mean 
field theory in this model. A discrepancy of this sort 
is a ubiquitous feature of mean field theories, which are 
principally designed to capture the mean field value. 

The comparisons so far have relied on the initial value 
A{t = 0) being of (negative) 0{K), that is, an initial dis- 
tribution whose width is of the same order as that of the 
steady state. When the initial width of the distribution 
is much larger than that of the steady state, i.e., when 
A[t = 0) is very different from (much smaller in magni- 
tude than) K, it becomes more problematic to capture 
the transient dynamics, although the steady state behav- 
ior is still predicted accurately. This is shown in Fig. [3l 

Finally and importantly, we note that the large cou- 
pling theory accurately captures the multistable nature 
of the dynamics (Fig. d]). This is a new feature of this 
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FIG. 2: (Color online) The time evolution of the mean 
(top panels) and inverse width (bottom panels) is shown 
for simulations of locally coupled two-dimensional arrays of 
TV = 64 X 64 sites and compared to the numerical solution 
of Eqs. (|20|l and (|2ip . The left column represents mod- 
est coupling, K — 10, while the right column shows re- 
sults for K = 80; cr^ = 3.5 for all plots. The dark (blue) 
lines show data from lattice simulations, while the light (red) 
lines represent theoretical predictions. The initial values are 
A{0) = —9.2 (left column) and -42 (right column). The inset 
in the bottom right panel shows a close up of the early time 
evolution. 
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FIG. 3: (Color online) The time evolution of the mean field 
(top panels) and the inverse width of the distribution (bottom 
panels) is shown for both locally (dark or blue) and globally 
(black) coupled simulations and compared to the numerical 
solution (light or red) of the large K theory. Again, the left 
column is for K — 10 and the right for K = 80. Here we 
choose 1^(0)1 < 1 (specifically, A{0) ~ —0.1), and the tran- 
sient dynamics are not well-described for early times. The 
theory accurately captures the steady-state behavior of the 
globally coupled simulations in the large K regime. = 4000 
(globally coupled) and 64 x 64 (locally coupled), and cr^ = 3.5 
for all panels. 
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FIG. 4: (Color online) Multistability is captured by the large 
K theory. The left panel shows the dynamics of the mean field 
and the inverse width of the distribution for initial conditions 
leading to a disordered state (iy£'o(0) « 0.5, A{Q) = -29). The 
right panel shows the corresponding plots for initial conditions 
leading to an ordered phase {(po{Q) ~ 4, A{Q) = —55). The 
inset in the bottom left shows a close up of the steady state 
behavior. Dark (blue): globally coupled simulations; light 
(brown): Gaussian ansatz theory; black: mean field theory. 
In all plots, K = 80, = 3.0. 



theory that provides dynamical information not provided 
by the usual mean field theory. In our earlier work we had 
established this multistability only through direct numer- 
ical simulations of the array. For large K, the agreement 
between theory and simulations both in the dynamical 
regime and in the steady state is quite remarkable for ini- 
tial conditions leading to the ordered state (right panel). 
While the mean field theory of course predicts the width 
of the distribution exactly, the Gaussian ansatz slightly 
overestimates the width for initial conditions leading to 
a disordered phase (left panel). The disordered state is 
marked by a relatively broad distribution, and, not sur- 
prisingly, given the underlying requirements stated ear- 
lier, the theory does not exactly capture this prediction. 



V. TWO FIELDS 

The more stringent test of the theory lies in the two- 
field system, where there are only numerical simulation 
results. Now we begin with Eqs. p?]) and p^ . substitute 
the assumed distribution into the two-variable Fokker- 
Planck equation (|13p . and again implement a large- if 
expansion. This results in the set of equations for the 
coefficients in the Gaussian llSll. 
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E = To{K + 2a'^A)E + 2uj{A - A/), 



(25) 



along with those for the mean values 



(26) 



Again, it is easy to ascertain that the normalization con- 
dition is consistent with the evolution equation ob- 
tained for C{t). 

The solution we are interested in is a collective oscilla- 
tory mode, which of course requires the time-dependent 
solution of the coupled sets of ordinary differential equa- 
tions ((25|) and ((26)) . a task which is vastly simpler than 
the solution of the time-dependent Fokker-Planck equa- 
tion. However, as in the single field case, we can further 
simplify the problem of finding the oscillatory long-time 
behavior of the mean values (/Sq and by exploring the 
regime where the coefficients A, M, and E have reached 
a steady state, that is, by setting the left hand sides in 
Eq. (|25p equal to zero. There are four stationary solu- 
tions, one of which is 



^, E = 0, M = A. 



(27) 



The other solutions 
(A, E, M) 
{A,E,M) = (0, 



{A, E, M) = (0, 0, 0), 
{-K/a^, -2w/(cr2ro), 0), and 
-2uj/{a'^To), -Kja^) do not sat- 



isfy the condition iMA — i? > necessary for proper 
normalization, cf. Eq. (|19p . and are hence unphysical. 

For the long-time behavior, it now remains to substi- 
tute the stationary A into Eq. ([^5]) and solve the coupled 
set of just two equations. 



dVeff((po) 

= -ro — , - t^zo, 



(28) 



the 



We note that this set has one stationary solution, 
disordered state (1^0,2:0) = (0,0). 

We can go even further toward the analytic oscillatory 
solution by implementing a multiscale perturbation the- 
ory For this purpose we combine l|26|) into a single 
second-order differential equation, 



z + uPz = ujG{lp), 



G(^) = -r(^) 



dip ' 



(29) 



with if = z/uj. It is understood that the variables are the 
mean fields and should therefore carry the subscript, 
which we have omitted for economy of notation. We treat 
G as a perturbation and write 



z + uP'z — eujG{z/Lo), 



(30) 



where e is a small parameter. The solution is then ex- 
pressed in terms of different time scales, (To, Ti, . . . ), 



z(t) = Z(To,Ti,...) 



Zo(ro,ri) + eZi(ro,ri) + o(62), 

(31) 
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where T„ 0(e"To). For the time derivatives we have 



dz 
'dt 



dz dz dz, 

dTo dTi dTo 



d^z 

9TS 



d^z 

'dTodTi 



(32) 



dTodTi 

The evolution equation at zero order in e is 



iiP'Zn 



0, 



which has a solution 

ZoiTo,Ti) = RiTi)cos{LuTo) + ^(Ti) sin(u;To). 



(33) 



(34) 



The coefficients R{Ti) and S{Ti) are determined by 
considering the evolution equations at first order in e, 



d^Zi 
on 



u?Zx = - 



dTodTi 



As always in multiscale perturbation theory, in order to 
avoid secular terms in the solution R and S must be 
chosen in such a way that resonant terms do not appear 
in the right hand side of Eq. ([55]) . i.e., 



dTo cos(a;To) 
dTo sin(a;ro) 



' dTodTi 
dTodTi 



-LuG 
ujG 



1 azo 

u dTo 
1 dZp 
uj dTo 



= 0, 



0, 

(36) 



where from ([M]) it follows that 
d'Zo OR 

Ll> sm(ci;Joj " 



dS 



ujcos{ujTo). (37) 



aTo^Ti dTi 
This immediately leads to the equations for R and 5, 

OR _ lZdToMu;To)G[i§P;) 



dTi 



/_oo sin (ujTo) 



(38) 



dS roo^roCos(c.To)G(lff) 
5Ti /_°°^dToCOs2(a;To) 

The quotients can be computed by introducing a cutoff 
A in the integrals, dTo ■ ■ ■ : E^nd subsequently taking 
the limit A — > oo. 

We thus have theoretical predictions at three levels of 
approximation. The most detailed is the five-equation set 
((25)) and (|26p . This set of equations contains the dynam- 
ical approach to the long time behavior. The second level 
is contained in Eqs. (P7)l and (I^S)) . This yields the long- 
time oscillatory behavior of the mean field and the width 



parameters of the distribution at long times. Finally, 
Eq. (H?]) together with ^ and ^ (along with the sec- 
ond equation in (|26p ) provide a full long-time solution 
in terms of quadrature. Note that this explicit solution 
presents a circular limit cycle of the form z — So sin(a-'t) 
with Sq constant [Zo ~ S'osin(a;ro) ], provided we have 
dR/dTi ~ [directly satisfied for the choice in and 
dS/dTi = 0, i.e., if 



= HSo) 



dTo cos{ujTa) G{So cos(cjTo)). (39) 



Since our numerical simulations [lyl indicate an essen- 
tially circular limit cycle near the onset of multistability, 
we expect that Eq. p9p may provide an accurate predic- 
tion of the onset in the strong coupling limit. 

We now proceed to test our multi-level theoretical pre- 
dictions against direct numerical simulations of the array. 
Again, we adopt the representative choices (|24p for the 
potential and relaxational functions. We concentrate on 
the multistability onset regime by fixing the value of the 
noise intensity appropriately. 

Figure \5\ compares the results of simulations for glob- 
ally coupled arrays with those of our theory. The agree- 
ment is clearly excellent for the limit cycle radius as well 
as the width parameters for all times. The time scale 
separation required for the validity of the large K equa- 
tions is clearly satisfied for K ^ 80 and even for modest 
values, K ~ 10. Note that the width parameters reach 
the steady state very quickly even when E is initially 
chosen to be of 0{K). The theory very slightly overes- 
timates the steady state radius r = \/ z? j (2w^) -|- z^/2 
for modest K = 10, and very slightly underestimates the 
radius for large X = 80 (neither visible on the scale of 
the figures). Again as expected, the agreement with the 
simulations for the locally coupled array is less spectac- 
ular (Fig. [6]) , but the limit cycle radius is still captured 
very accurately for all times. For visualization purposes, 
we also show the early time evolution of the theoretical 
probability distribution function in Fig. [71 

Figure [5] confirms that the phase portrait of the limit 
cycle from the globally coupled lattice simulations cor- 
responds reasonably well to that predicted by the large 
K equations, and that the frequency of the oscillations 
is also accurately predicted. The dynamic evolutions of 
the simulation and of the theory do not match exactly 
because the small error in the frequency implies an in- 
creasing dephasing with time. The figure shows the case 
K = 80, but the results are only representative, and 
again the theory holds well even for more modest val- 
ues of K and for a range of initial conditions leading to 
limit cycle behavior. 

Finally, as in the single field problem, the large cou- 
pling theory predicts the occurrence of multistable re- 
gions in parameter space. In our earlier work [l0| we 
noted the expectation of multistability, but did not ex- 
plicitly pursue it in our simulations, nor did we have a 
predictive theory as we now do. In Fig. [5] we show that 
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FIG. 5: (Color online) The time evolution of the limit cycle 
radius, r (top panels) and of the Gaussian ansatz coefficients 
(three subsequent panels) . The dark (blue) curves are appro- 
priate moment results from simulations of globally coupled 
arrays {N — 4000 sites), and the light (red) curves are ob- 
tained from the solution of Eqs. (|25p and (|26p . = 10 in the 
left column and if = 80 in the right column. Initial values 
for left column: {A,E,M) = (-9.6,4.7,-9.7); right column; 
{A, E, M) = (-56, 7, -79). The noise intensity = 3.5 in aU 
panels. The dashed line in the top panels is the steady state 
radius predicted by multiscale analysis. 
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FIG. 6: (Color online) Same as Fig. O but now the sim- 
ulations are for locally coupled arrays. Initial values for 
left column: {A,E,M) = (-9.4,4.4,-9.6); right column: 
{A,E,M) = (-55,5,-76). 



FIG. 7: (Color online) Evolution of the theoretical probability 
distribution function as obtained by substituting the numer- 
ical solution of Eqs. ([25} and JMl into Eq. The distri- 
bution becomes symmetrical (circular) in i^o and zq prior to 
reaching its steady state radius. The evolution is only shown 
for the short time leading up to relaxation of the distribution 
shape. Over longer times, the distribution will continue a 
circular trajectory whose radius eventually reaches its steady 
state value (see Fig. [8]). K = 10 and — 3.5. 

results indicate that at if = 10 multistability first oc- 
curs at roughly cr^ 3, with the noise value decreasing 
ever more slowly with increasing K (e.g., cr^ « 2.6 for 
K = 20, « 2.5 for K = 40). At if = 80 the transition 
occurs at « 2.4. Equation is in fact independent 
of K. The function I{So) is shown explicitly in Fig. [TOl 
The first zero is the disordered state and the second first 
appears when <t^ = 2.4. As noise increases, an intermedi- 
ate unstable solution also appears, illustrated explicitly 
in Fig. [9l The multiscaling result thus accurately predicts 
not only the radius r = 5o of the limit cycle but the noise 
for onset of multistability when coupling is strong, and 
the increase in the limit cycle radius with increasing noise 
strength. 



the multistable nature of the phase transition is fully cap- 
tured by the large coupling theory. In particular, differ- 
ent initial conditions lead to either a disordered phase 
or a limit cycle. To test whether the theory correctly 
predicts the onset of multistability is more cumbersome, 
but we can at least do it easily in the multiscaling (large 
K) regime on the basis of Eq. p9|) . The value of the 
noise for onset of multistability in general depends on 
coupling strength [l3|, but this dependence greatly 
weakens with increasing K. For example, our simulation 



VI. CONCLUSIONS 

We have presented an analytical theory for the dy- 
namics of relaxational systems with field dependent co- 
efficients. Based on a Gaussian ansatz and an expansion 
about the mean field values, we derive ordinary differ- 
ential equations detailing the time evolution of the field 
distribution means and widths. In the limit of large cou- 
pling, the Gaussian ansatz equations provide a consis- 
tent, normalized approximation of the relevant probabil- 
ity distributions which agree with numerical lattice simu- 
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FIG. 8: (Color online) The phase space portraits of lattice 
simulations (dark or blue) and the large K equations (|25|l and 
(|26|1 (light or red) show good agreement in terms of the steady 
state limit cycle radius, though the theory slights underesti- 
mates r. Inset: oscillation frequency. K = 80, — 3.5, and 
= 1. 




0.00.1 



lations. In particular, our method allows us to study the 
dynamics of both one- and two-field relaxational systems, 
including those with oscillatory collective states previ- 
ously beyond the reach of the static mean field theories. 
This also provides an analytic access to initial-condition- 
dependent multistable regimes for noise-induced phase 
transitio ns in spatially extended systems. 




FIG. 10: The multiscale analysis predicts the onset of bista- 
bility to occur at ~ 2.4. Curves are, from thickest to 
thinnest, cr^ = 3.5, 3.0, 2.4, 2.0. 
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FIG. 9: (Color online) Comparison of the phase space por- 
traits of lattice simulations (dark or blue) and the large K 
equations (light or red) show that the latter correctly predict 
multistable behavior. The right inset shows a close up of the 
indicated portion of the phase portrait. The black curve is 
the unstable limit cycle obtained from the multiscale analy- 
sis. K = 80, — 3.5, and uj — 1. 
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